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摘 要 建立 了 宽 弦 叶片 定向 凝固 过 程 的 宏 微 观 数理 模型 , 计算 结果 与 实验 结果 吻合 , 铸件 表面 大 部 分 晶 粒 的 位 置 和 形状 

对 应 . 通过 数值 模拟 预测 了 不 同 引 晶 方 式 和 拉 速 下 温度 场 和 唱 粒 组 织 的 演变 过 程 , 研究 了 2 种 因素 的 影响 规律 . 通过 建 
开 糊 状 区 形态 和 唱 粒 数 的 数学 判 据 , 实现 了 温度 场 和 微观 组 织 优 务 的 定量 评价 , 基于 判 据 揭示 了 工艺 参数 对 糊 状 区 和 晶 粒 
的 影响 机 理 , 从 而 对 工艺 进行 了 定量 的 优化 . 研究 表明 , 采用 叶 身 正 下 方 的 引 晶 方式 , 有 利于 增加 柱状 晶 数 目 , 提高 晶 粒 平 
行 度 , 防止 横向 晶 界 生成 , 同时 还 可 以 在 糊 状 区 形状 保持 平 直 的 前 提 下 使 用 较 高 拉 速 , 从 而 避免 晶 粒 粗 化 , 并 提高 生产 率 . 
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ABSTRACT The aero turbine is spun by high-temperature and high-pressure burning gases. The practice has 
proven that the directional solidification (DS) turbine blade with perfect column grains has still excellent high-tem- 
perature performance in this kind of working environment. This means that the size and orientation of column 
grains have great influence on the high-temperature property and performance of turbine blades. On the other hand, 
the high-quality blade is not easy to be produced in DS process due to the difficulty of obtaining the desired temper- 
ature field needed to produce the grains with ideal morphology. In addition, the growth of columnar grains in the 
wide- chord hollow guide blade is obstructed by the complex camber and the platform. How to produce turbine 
blades with desired microstructures is the key problem in the DS process. Numerical simulation of the DS process 
is an effective way to investigate the growth and the morphology of the grains and hence to optimize the process. 


In this work, a mathematical-physical model for simulating the DS process of wide-chord blade is established in 
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which nucleation and grain growth in the blade in the DS process are modeled by the cellular automation (CA) 
method with multi-scale dynamic bidirectional coupling technology. Some general analytic indicators are proposed 
to assess the morphology of mushy zone and grains in a blade quantitatively. Based on the simulated results by us- 
ing the usual starter blocks 1, 2 and 3, a new starter block is designed considering numerically controlled cutting. 
Temperature fields and grains in DS processes and corresponding indicators at different withdrawal rates for above 
4 starter blocks are numerically predicted to investigate the influences of varying these technological parameters, 
and hence to determine the influence mechanism to the DS process. For comparison, the DS validation experiments 
by using starter blocks 1, 2 and 3 have been carried out. The numerical and experimental results agree well, their 
morphologies including those faulty grains are similar. It is found that higher withdrawal rate leads to larger conca- 
vation of mushy zone, but the effect of chill is stronger than that of withdrawal rate if the contact area between cast- 
ing and chill plate is large enough. Better grain structure in a blade is achieved by starter block 3 than by starter 
blocks 1 and 2. By starter block 4, the amount of column grains is larger and the amount of lateral grain boundaries 
is smaller, compared with that of starter blocks 1, 2 and 3. Therefore higher withdrawal rate could be adoptable 
without excessive concavation of mushy zone, resulting in parallel column grains, finer dendrites in the blade, and 
much higher blade productivity. Optimum withdrawal rates are also determined for starter blocks 3 and 4. 


KEY WORDS aero turbine, wide-chord blade, directorial solidification, numerical simulation, columnar grain 


我 国航 空 工业 的 发 展 急需 要 提高 航空 发 动机 ” 固 过 程 一 般 采 用 Bridgman 法 上 ", 要 求 对 温度 场 
的 综合 性 能 ". 高 压 涡轮 导向 叶片 的 高 温 力 学 性 能 。 行 极 划 ed 包括 拉 速 、3 
和 持久 性 能 是 发 动机 性 能 的 决定 因素 之 一 , 叶片 的 唱 方 式 等 , 定向 凝固 工艺 窗口 很 难 寻找 . 用 实验 i 
柱状 唱 ， 其 苛刻 的 限定 中 例如 , 偏离 度 行 工艺 研 究 的 周期 长 .成 本 高 . 因此 , 很 有 必要 应 用 


°, 发 散 度 应 不 大 于 25?; 叶 身 中 间 柱 状 。 ”数值 模拟 技术 来 对 温度 场 进行 辅助 研究 . 
少 、 蝇 粒 不 能 大 宽 . 这 些 都 是 为 了 减 近年 来 , 人 们 开始 采用 计算 机 数值 模拟 技术 研 
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少 横向 晶 界 , 特别 是 叶 身 进 、 排 气 边缘 上 不 允许 存 。 究 高 温 合金 定向 凝固 微观 组 织 W Yang 等"… "和 
在 横向 唱 界 , 此 外 , 还 要 求 微观 组 织 细密 , 具有 较 小 。” Wang 等 "用 元 胞 自动 机 (CA) 法 模拟 了 高 温 合金 定 
的 校 唱 臂 间距 . 但 是 , 高 压 导 向 叶片 具有 和 宽 弦 、 大 缘 。” ”向 凝固 试 件 的 局 部 校 蝇 演变 , 并 以 此 解释 杂品 形成 
板 、 注 壁 和 复合 倾斜 的 形状 特点 , 如 图 1 所 示 . 这 些 。 机 理 . Gandin 等 只 和 Rappaz 等 创建 了 基于 枝 唱 
结构 给 叶片 定向 凝固 过 程 中 的 柱 晶 生长 带 来 很 大 尖 包 络 推 移 的 3D-CA 算法 , 使 得 计算 较 大 区 域 的 品 


的 干扰 . 晶 粒 从 下 缘 板 引 品 段 生 长 后 转向 , 罕 过 横 


粒 形 貌 成 为 可 能 . Dai 等 中 和 Yang 等 外 基于 商用 软 


向 缘 板 到 叶 身 , 再 次 转向 进入 叶 身 . 截面 的 变化 容 。 件 ProCAST 的 元 胞 自动 机 -有 限 元 (CAFE) 模 块 模拟 
易 引 起 晶 粒 粗大 、 断 晶 以 及 晶 粒 取向 偏离 度 大 等 缺 了 单 唱 叶片 在 Bridegman 下 的 引 唱 器 、 选 晶 器 唱 粒 组 
陷 ? 53. 另 一 方面 , 目前 高 压 涡轮 导向 叶片 的 定向 凝 ” 织 . Gandin 等 模拟 了 航空 定向 凝固 叶片 的 晶 粒 组 


图 1 宽 弦 导 
Fig.1 3D geometry of wide-chord guide blade 


向 叶片 的 三 维 结构 


织 , 但 未 进行 实验 验证 .Kermanpur 等 所 用 商用 软件 
ProCAST 的 CAFE 模块 计算 了 燃气 轮机 定向 凝固 中 


工 


2 晶 粒 组 织 , 预测 了 特定 工艺 下 的 唱 粒 发 散 生 


长 , 但 没有 准 


确 预 测 到 实验 中 的 断 晶 和 晶 


粒 宽度 分 


布 . Pan 等 82 
了 定向 凝固 


Bridgman 航空 


基于 自主 
的 网 格 量 的 问题 , 进而 在 国 


发 的 软件 , 用 分 层 算法 解决 


内 率先 对 


涡轮 叶片 全 铸件 晶 粒 组 


织 进行 了 计 


算 ; 尚 伟 等 的 利用 ProCAST 计算 了 单 晶 导向 叶片 的 


晶 粒 组 织 , 预测 了 杂 晶 ; ee p 用 ProCAST 计算 
了 单 晶 试 板 的 微观 组 织 ,在 


究 了 拉 速 和 缘 板 大 小 对 


糊 状 区 形状 的 影响 , 但 没有 研究 糊 状 区 形状 随 高 度 


( 抽 拉 距离 ) 的 变化 , 因 


晶 的 影 响 . 


叶片 柱状 草 或 单 唱 组 织 的 预测 , 需 


而 无 法 衡量 缘 板 的 高 度 对 杂 


需要 通过 求解 
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全 铸件 的 微观 网 格 来 实现 . 分 层 计 算法 ”在 每 个 时 
间 步 只 对 糊 状 区 及 其 附近 的 微观 组 织 进行 求解 , 从 
而 节省 内 存 和 CPU 时 间 , 但 是 , 在 定向 凝固 过 程 中 ， 
叶片 的 宽度 决定 了 凝固 前 沿 的 面积 , 宽 弦 叶片 的 糊 
状 区 体积 较 大 , 因此 , 分 层 计 算 在 每 个 时 间 步 消耗 的 
内 存 和 时 间 也 就 大 , 叶片 的 宽度 成 为 计算 的 瓶颈 . 

目前 , 对 叶片 定向 凝固 晶 粒 组 织 的 模拟 存在 以 
下 问题 : (1) 计算 尺度 太 小 , 计算 的 试 板 和 叶片 主体 
部 分 宽度 都 在 50 mm 以 内 ; (2) 对 试 板 的 研究 较 多 ， 
也 有 对 涡轮 叶片 的 模拟 , 但 是 对 宽 弦 导向 叶片 的 模 
拟 很 罕见 ; (3) 对 计算 结果 的 可 视 化 观察 过 程 缺乏 定 
量 的 判 据 , 特别 是 当 不 同 工 艺 下 的 计算 结果 大 同 小 
异 时 , 无 法 对 组 织 优 劣 进行 更 精确 的 比较 ; (4) 对 晶 
粒 组 织 形 貌 的 计算 , 从 总 体 趋 势 和 统计 规律 上 已 经 
能 与 实验 吻合 , 但 在 随机 形 核 模型 影响 下 , 具体 每 
个 晶 粒 的 位 置 和 形状 很 难 与 实验 一 一 对 应 , 尤其 是 
较 受 关注 的 特殊 晶 粒 , 如 断 晶 、 分 又 晶 等 ; (5) 国内 的 
研究 依赖 于 ProCAST 等 国外 商用 软件 , 一 些 特殊 的 
计算 需求 (如 柱状 晶 定量 比较 判 据 ) 无 法 实现 . 
本 工作 基于 自主 开发 的 铸造 软件 TecCast, 对 定 
向 凝固 过 程 中 温度 场 和 晶 粒 演变 进行 建 模 和 动态 
双向 耦合 计算 , 并 以 此 实现 宽 纺 航空 叶片 的 全 铸件 
晶 粒 计算 , 进而 了 解 其 演变 机 理 ; 建立 相关 数学 判 
据 , 对 不 同 工 艺 因素 下 糊 状 区 形态 、 唱 粒 发 散 度 进 
行 定量 评价 和 优化 . 如 前 所 述 , 宽 弦 导 叶 晶 粒 组 织 
不 合格 的 主要 原因 是 , 晶 粒 从 引 晶 段 生 长 到 叶 身 
时 , 未 能 形成 优化 的 柱状 晶 粒 . 另外 , 在 定向 凝固 
中 , 拉 速 是 核心 参数 之 一 , 因此 , 本 工作 从 引 晶 方式 
和 拉 速 2 个 方面 进行 数值 计算 和 工艺 优化 . 
1 宽 弦 叶片 定向 凝固 过 程 宏 微观 数学 模型 
图 2 可见 , 高 温 合 金 浇 注 后 , 型 壳 随 水 冷 Cu 
板 一 起 自 上 而 下 运动 . 在 定向 凝固 过 程 中 , 辐射 是 
最 主要 的 热 交换 方式 . 加 热 区 炉 壁 对 叶片 起 辐射 加 
热 作用 . 当 模 组 进入 冷却 区 后 , 叶片 模 组 与 冷却 区 
水 冷 炉 壁 间 发 生 剧 烈 的 辐射 散热 . 在 这 个 过 程 中 ， 
在 叶片 中 形成 自 上 而 下 的 温度 梯度 , 完成 定向 凝 
固 . 本 工作 基于 VC++ 编 制 软件 , 对 温度 场 和 唱 粒 形 
核 生 长 的 数学 模型 进行 求解 计算 . 
1.1 温度 场 

温度 场 控 制 方程 : 

六 


oT _ ,mo of, 
pov dr -AV TtpAH ES + Ora (1) 


式 中 , p 为 密度 , kg/mi; c 为 比热容 , Jkg.K); 7 为 热 
力学 温度 , K; 1 为 时 间 , s; 4 为 导热 系数 , W/(mK); 
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Heating 


e 
Chill plate 


和 
SS 


2 Bridgman 法 示意 图 


Fig.2 Schematic of Bridgman method 


A 为 结晶 潜 热 , J/kg; 为 固 相 率 ; O, 为 实体 表面 单 
元 与 环境 之 间 的 热 交 换 , W/m:. 
实体 表面 辐射 传 热 采用 改进 的 Monte Carlo 射 
线 追 踪 法 "来 计算 . 温度 场 边 界 条 件 为 : 
BT 于 7 


_ n 2 
Qi > ll —e (1 e,)S, | | 
uw +1 


&, £5, 

式 中 , Oo 为 型 壳 或 结晶 盘 表 面 单元 的 自由 表面 
与 环境 之 间 的 热流 密度 , W/m; $, 和 5, 分 别 为 u 的 面 
积 和 射线 vy 的 环境 子 空间 的 面积 , m2; o 为 Stefan- 
Boltzmann 常数 , W/m2.K9; Bi, 为 改进 的 射线 y 的 权 
重 系数 ; n 为 u 发 出 的 射线 总 数 ; 雹 和 工分 别 为 uw 和 vw 
的 温度 , KK; &, 和 sg, 分 别 为 w 和 v 的 灰 度 . 其 中 v 有 可 能 
是 热 区 和 固定 挡 板 表面 , 也 可 能 是 型 这 或 结晶 盘 的 
另 一 个 自由 表面 . 每 一 个 单元 的 空间 坐标 ZZ, 在 每 一 
个 时 间 步 根据 拉 速 实时 更 新 , 每 隔 10 个 时 间 步 更 
新 7. 
1.2 微观 组 织 

1.2.1 形 核 模型 采用 Rappaz 和 Gandinc 的 
Gauss 分布 随机 形 核 模型 : 

oN  _- NN on- (A7,- 2 (3) 
AAT) VanAT, 2(AT,) 

式 中 , 和 N 为 形 核 密度 , m’; AT 为 当前 热力 学 过 冷 度 ， 
K; N 为 最 大 形 核 密度 , mi; A7, 为 形 核 过 冷 度 标准 
方差 , K; AT 为 平均 形 核 过 冷 度 , K. 其 中 MX AT 和 
ATx 结 合 引 晶 块 的 唱 粒 观察 结果 取经 验 值 . 

于 是 , 晶 粒 密度 MAT) 可 表示 为 : 


才 
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NAT)=h FT) (0) 
式 中 ，A7" 为 过 冷 度 积 分 元 . 
1.2.2 唱 粒 生长 模型 
取 近 似 的 KGT™"" 模 型 , 即 : 
v=aAT? +BAT’ (5) 
式 中 , v 为 枝 晶 尖 推 移 速 率 , m/s; w 和 为 生长 系数 ; 
A7 为 枝 晶 尖 生 长 前 沿 的 过 冷 度 , K. 
凝固 过 程 中 校 晶 生长 前 沿 的 过 冷 度 由 4 部 分 
组 成 : 


曲 粒 生长 模型 可 以 采 


AT=AT+AT.+AT +AT (0) 

式 中 ，A7. 为 成 分 过 冷 度 ，AT7, 为 曲率 过 冷 度 ，AT 
为 动力 学 过 冷 度 . 其 中 动力 学 过 冷 度 一 般 远 小 于 前 
3 项 , 可 以 忽略 不 计 . 
1.3 宏 微观 动态 双向 耦合 

采用 双重 嵌 套 网 格 , 宏观 网 格 用 于 温度 场 计 
算 , 每 个 宏观 网 格 中 划分 为 若干 微观 网 格 , 用 于 微 
观 组 织 的 CA 捕获 计算 . 在 以 往 CAFE 全 叶片 唱 粒 计 
算 中 , 为 节省 计算 量 , 研究 者 大 都 采用 宏 微 观 间接 
耦合 , 亦 即 仅 从 宏观 温度 场 插值 到 微观 温度 场 的 单 
向 灯 合 “为 了 更 加 精确 地 计算 枝 晶 潜 热 释放 对 
温度 场 的 影响 , 本 工作 采用 双向 动态 同步 耦合 . 

温度 场 的 一 个 时 间 步 完成 后 , 将 温度 值 插值 到 
微观 网 格 : 


Pe = T+ Cr “6x— %) 


Gr(m*6y— 5 )+Gr(n:6e— %) (7) 
式 中 ,i,j 和 分 别 为 宏观 单元 在 x,y 和 z 方 向 上 的 序 
号 ; 1,m 和 nn 分 别 为 宏观 单元 内 部 坐标 系 下 微观 网 格 
在 x,y 和 z 方 向 上 的 序号 ; ee 为 微观 网 格 的 温度 ， 
K; 人 六 为 宏观 单元 的 温度 , K; G, G; 和 Gi 分 别 为 x,y 
Hz 方向 上 的 温度 梯度 , K/m; bx, 6y 和 6z 分 别 为 x,y 
1z 方 向 上 的 微观 网 格 尺寸 , m; a, b 和 cc 分别 为 x,y 
1z 方 同上 的 宏观 单元 尺寸 , m. 
当 1.6x 二 04, G 取 向 前 差分 , 反之 向 后 差分 , G 
[Gi 同 理 . 
该 时 间 步 对 应 的 微观 计算 完成 后 , 将 微观 节点 
的 烩 变 累 计 到 宏观 节点 : 
A 厅 入 =D， AH'Y, (AR Bx “6y*6z) (8) 


A 


污 


E23 


个 宏观 时 间 步 的 固 液 


式 中 , A 有 ,,, 为 宏观 单元 在 


相 变 潜 热 ,J，AF4"" 为 微观 网 格 在 宏观 时 间 步 的 固 
相 率 变 化 . 


ChinaXiv 合 作 期 刊 


学 报 第 51 郑 


2 凝固 前 沿 形 瑶 与 晶 粒 数 的 数学 判 据 

为 了 定量 衡量 温度 场 对 微观 组 织 的 影响 以 及 
不 同 工 艺 下 微观 组 织 的 优 劣 , 对 于 不 同 工 艺 下 得 到 
的 糊 状 区 形状 和 柱状 晶 生 长 状况 , 需要 一 个 数值 化 
的 标准 来 评判 , 以 定量 衡量 工艺 的 优 劣 , 为 其 优化 
提供 精确 的 客观 依据 . 

2.1 糊 状 区 宽度 

对 于 准 稳 态 定向 凝固 , 一 次 枝 唱 臂 间距 随 凝 固 
前 沿 温度 梯度 递增 , 该 梯度 近似 与 糊 状 区 宽度 成 反 
比 ; 二 次 枝 晶 臂 间距 随 局 域 固 液 共存 时 间 递 减 , 该 
时 间 近 似 与 糊 状 区 宽度 成 正比 . 因此 , 糊 状 区 宽度 
越 小 , 枝 唱 越 细密 . 

某 一 时 刻 的 糊 状 区 平均 宽度 到 (0 为 : 

及 (D=T/S (9) 
式 中 , ,为 等 温 面 和 7 等温 面 之 间 包 围 的 铸件 
部 分 的 体积 , m’; 和 和 也 分 别 为 铸件 的 固 相 线 温度 
和 液 相 线 温度 , K; $ 为 该 体 在 xoy 平 面 上 的 投影 面 
避 , mm2， 
2.2 凝固 前 沿 下 凹 度 

由 枝 晶 生长 模型 可 知 , 梳 晶 尖 处 的 温度 梯度 的 
方向 是 唱 粒 生长 方向 的 主要 影响 因素 , 因此 , 糊 状 
区 凹凸 决定 柱状 晶 发 散 度 , 如 图 3 所 示 . 一 般 情况 
下 史 , 对 于 一 定 尺寸 的 铸件 , 抽 拉 太 慢 会 导致 糊 状 区 
凸 度 增 大 , 柱状 晶 产 生 V 形 偏向 , 亦 即 向 上 发 散 ; 抽 
拉 太 快 会 导致 糊 状 区 凹 度 增 大 , 柱状 唱 产 生 A 形 偏 
向 , 亦 即 向 下 发 散 ( 反 发 散 ) 叫 . 通过 优化 拉 速 可 以 避 
免 糊 状 区 凹凸 , 提高 晶 粒 的 平行 度 . 

一 般 定 向 凝固 的 情况 下 , 枝 晶 尖 离 也 等 温 面 都 
很 近 (1 mm 以 内 ), 因此 7 等 温 面 可 近似 看 作 糊 状 区 
包 络 面 , 计算 中 近似 用 7 等温 面 上 各 点 的 温度 梯度 
来 表征 各 枝 晶 尖端 的 温度 梯度 . 

在 图 4 中 , 某 时 刻 铸件 合金 部 分 的 等温 面 为 
了 界 曲面 2 显示 为 红色 , 其 边界 在 型 过 人 
上 , 其 在 xoy 平 面 上 的 投影 为 8' 投影 面积 为 9, m 
0 为 该 时 刻 8 的 几何 中 心 点 ; Zo 为 竖 直 向 上 的 无 量 
纲 单 位 矢量 ; P 为 0 上 任 一 点 , 其 邻 域 为 dx dy, m; 
Gs 为 点 P 处 的 温度 梯度 , K/m, 其 方向 沿 8 在 点 P 处 
的 法 癌 . 

点 O 的 定义 满足 : 

for=0 (10) 


该 时 刻 糊 状 区 下 凹 度 K。.(D 为 : 


0=] (OP x2,)(G, x 2,) 
» SIoPlG,| 


dxdy (11) 


2 
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(a) 


Liquid 
Liquid 
Liquid 


图 3 拉 速 对 凝固 前 沿 轮廓 的 影响 


Fig.3 Profile of solidification front dependent on the withdrawal rate (G—temperature gradient at solidification front) 


(a) excessively fast (b)moderate rate (c)excessively slow 


(b) 


Zo 


图 4 凝固 前 沿 下 凹 度 的 计算 方法 示意 图 
Fig.4 3D (a) and 2D (b) schematic diagrams of evaluation methodology used for calculating the concavation of solidifica- 
tion front in Eq.(11) (2 一 isothermal surface of 7T'; O 一 geometric centre point of 2; Zo 一 vertical up unit vector; P— 


point on 2 with neighborhood dxdy; G,—temperature gradient along the normal direction of P) 


Kw>0 显示 联 等 温 面 下 凹 , Ku<0 显示 了 等温 薄 壁 叶片 , 表面 横 截 线 晶 粒 数 能 间接 反映 横 截面 晶 
上 凸 .优化 工艺 的 主要 目的 之 一 是 减 小 K. 的 绝对 值 。 。” 粒 数 . 

2.3 横 截面 晶 粒 数 在 图 5 中 , 对 于 高 度 为 z 的 水 平 截面 与 铸件 表 
横 截面 晶 粒 数 的 变化 可 在 一 定 程 度 上 反映 柱 ” 面 形 成 横 截 线 L, 因 工 为 封闭 曲线 , 表面 模 截 线 蝇 粒 
状 唱 竞争 生长 及 生长 方向 的 变化 . 在 图 5 中, 与 高 度数 表示 为 了 与 唱 界 的 交点 数 , 定义 为 Ksvo(z"). 


为 = 的 水 平 截面 相交 的 晶 粒 数量 , 亦 即 横 截面 上 可 目前 商用 软件 中 必 ee 的 计算 程 
见 的 唱 粒 截面 的 个 数 , 定义 为 Ko(2). 序 , 本 工作 基于 VC++ 编 制 后 处 理 模块 , 对 温度 场 和 


2 4 表面 模 截 线 晶 粒 数 唱 粒 计算 结果 进行 数据 处 理 , 给 出 判 据 值 . 


在 实际 生产 中 , 检查 表面 晶 粒 是 比 检查 横 截 面 。 3 实验 验证 
晶 粒 更 为 快速 的 方法 , 已 有 对 应 的 工业 标准 对 于 铸件 采用 新 型 金属 化 合 物 型 1C10 合金 ". 采用 
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2 叶片 模 组 , 轴 对 称 分 布 , 如 图 2 所 示 , 壳 厚 7 mm. 初 
始 时 结晶 盘 上 表面 与 挡 板 上 表面 高 度 一 致 , 型 壳 预 
热 温 度 和 浇注 温度 均 为 1520 'C, 浇注 后 静 置 3 min 
开始 抽 拉 , 抽 拉 过 程 中 定向 炉 热 区 加 热 温 度 为 
1500 'C. 由 于 采用 顶 注 , 绝 大 部 分 金属 液 不 会 流 经 
结晶 盘 , 浇注 过 程 中 金属 液 温度 变化 不 大 ; 对 类 似 
过 程 的 计算 “表明 , 静 置 过 程 中 金属 液 温 度 趋 于 
均匀 . 

为 了 验证 计算 结果 , 对 不 同 工 艺 进 行 了 实验 , 采用 
90%HCI+10%H2O0;( 体 积分 数 ) 溶 液 腐蚀 晶 粒 组 织 . 
4 结果 与 讨论 
4.1 基于 数值 模拟 的 引 晶 方式 优化 

为 了 克服 下 缘 板 对 晶 粒 的 影响 , 必须 采用 引 唱 


图 5 晶 粒 数 的 计算 方法 示意 图 


Fig.S Schematic diagram of evaluation methodology used 


for calculating grain number (ZL 一 intersection line 


of casting surface and plane at z") 


结构 以 使 晶 粒 进入 叶 身 . 原 有 的 工艺 有 3 种 引 晶 方 
式 : 引 晶 1, 仅 在 下 缘 板 凸 楼 处 引 蝇 , 利用 了 零件 本 
身 设 计 的 向 下 延伸 结构 ; 引 唱 2, 在 引 晶 1 的 基础 上 ， 
在 下 缘 板 边缘 一 侧 添加 引 唱 块 ; 引 唱 3, 在 引 晶 1 的 
基础 上 , 在 叶 身 下 端面 , 靠 叶 盆 弱 线段 中 间 位 置 加 
引 唱 块 . 

引 晶 1,2 和 3 如 图 6a~c 所 示 . 根据 常规 设计 原 
则 , 引 唱 1 2 和 3 便于 成 形 后 切割 引 唱 块 . 在 拉 速 为 
7 mm/min 的 条 件 下 对 3 种 引 晶 方式 进行 了 数值 模 
拟 , 计算 得 到 各 个 时 间 步 的 温度 场 和 唱 粒 组 织 , 基 
于 计算 结果 , 用 前 文 所 述 的 算法 得 到 各 判 据 值 , 以 
考察 不 同 引 唱 方 式 对 温度 场 和 晶 粒 的 影响 . 对 3 种 
不 同 的 引 晶 方式 进行 了 实验 验证 . 随后 , 通过 分 析 
引 唱 方式 的 影响 规律 和 原 有 引 晶 方式 的 缺点 , 提出 
了 新 的 引 晶 方式 4, 如 图 6d 所 示 . 引 晶 4 为 在 引 晶 1 
的 基础 上 , 以 叶 身 下 端面 整体 向 下 延伸 至 结晶 器 . 
此 方案 需要 设计 专门 的 模具 , 并 采用 数控 铣床 或 线 
切割 床 进行 后 期 切割 . 对 引 唱 4 在 7 mm/min 工 艺 下 
的 温度 场 和 晶 粒 进行 了 模拟 . 

4.1.1 不 同 引 晶 方式 的 糊 状 区 形态 分 析 ”图 7 
为 不 同 引 晶 方式 叶 身 KK 的 变化 . 图 中 横 坐 标 为 凝固 
前 沿 到 激 冷 面 的 距离 , 叶 身 部 分 对 应 40~100 mm 高 
度 段 . 可 见 , 随 着 凝固 前 沿 远 离 结 晶 器 , K, 逐 渐 增 
大 . 这 是 由 于 激 冷 盘 单 位 面积 上 的 冷却 效率 较 高 ， 
而 模 壳 对 冷 区 辐射 换 热 时 , 单位 面积 上 的 冷却 效率 
比 激 冷 盘 低 得 多 , 因此 随 着 糊 状 区 到 激 冷 盘 之 间 的 
距离 增 大 , 糊 状 区 中 竖 直 方向 上 的 温度 梯度 下 降 . 


图 6 4 种 不 同 的 引 晶 方式 


Fig.6 Four different configurations of starter blocks 


(a)starter ]  (b) starter 2 


(c) starter 3  (d) starter 4 


EI 
SU 


202303.00464v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


第 4 期 再 宁 等 : 宽 弦 航空 叶片 Bridgeman 定向 凝固 组 织 数 值 模拟 505 


在 接近 上 缘 板 时 , K 有 一 定 减 小 , 这 是 因为 在 糊 状 
区 下 方 的 冷却 速率 保持 连续 的 同时 , 上 方 模 截 面 的 
忽然 增 大 , 使 得 凝固 前 沿 推 移 所 需要 的 热流 增 大 ， 
准 稳 态 热流 被 打破 , 糊 状 区 上 方 的 冷却 速率 突然 减 
小 , 使 糊 状 区 内 竖 直 方向 上 的 温度 梯度 突然 上 升 . 
相对 于 引 晶 1, 引 晶 2 和 引 唱 3 的 合金 与 激 冷 盘 接触 
面积 有 所 增加 , 提高 了 冷却 效率 , 因此 其 瑟 比 引 唱 1 
的 小 . 引 晶 2 与 引 唱 3 相 比 , 引 晶 2 的 引 晶 块 和 叶 身 
之 间 经 过 横向 且 较 薄 的 缘 板 连接 , 而 引 晶 3 的 引 唱 
块 直接 连接 叶 身 , 因而 引 晶 3 的 冷却 效率 更 高 , Ks 更 
小 .但 是 整体 上 引 晶 1,2 和 3 的 KK 差别 并 不 大 . 引 昂 
4 大 大 增加 了 与 激 冷 盘 的 接触 面积 , 且 直 接连 接 叶 
身 , 因此 , 其 Ki 明显 比 前 3 种 引 唱 小 , 这 有 利于 得 到 
细密 的 枝 唱 . 

图 8 为 不 同 引 唱 方 式 叶 身 的 及 变化 . 总 体 来 
说 ,4 种 工艺 的 凝固 前 沿 都 经 历 了 先 上 凸 再 下 目的 
过 程 . 在 铸件 靠近 模 壳 的 部 分 (铸件 表层 ), 热流 往往 
受 模 充 辐射 的 影响 较 大 , 其 凝固 前 治 的 位 置 高 度 主 
要 取决 于 挡 板 位 置 , 因而 激 冷 面 的 传 热 对 铸件 心 部 


10 
g| 一 Starter 1,7 mmM/min 
---- Starter 2, 7 mm/min 
CA Starter 3, 7 mm/min 


一 Starter 4, 7 mm/min 


本 
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Fig.7 Variations of thickness of mushy zone Ku in blade 


body with height for different starters 


0.4 
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8 不 同 引 晶 方 式 下 叶 身 糊 状 区 下 四 度 K 沿 高 度 的 

变化 


Fig.8 Variations of concavation of solidification front Ke 


in blade body with height for different starters 


凝固 前 沿 高 度 的 影响 大 于 对 铸件 表层 的 影响 . 在 叶 
身 下 半 段 , 亦 即 凝固 前 沿 离 激 冷 面 较 近 时 , 铸件 心 
部 的 凝固 前 沿 推移 较 快 , 容易 超过 表面 凝固 前 沿 
的 高 度 , 造成 整个 凝固 前 沿 上 凸 . 反之 , 叶 身 上 半 
段 凝固 前 沿 下 四 . 对 于 叶 身 底 端 , 由 于 激 冷 面 使 得 
模 壳 的 底 端 温度 较 低 , 铸件 表层 凝固 前 沿 高 度 不 
完全 取决 于 挡 板 位 置 , 因而 底 端 的 上 凸 程度 较 小 . 
对 于 叶 身 顶端 , 由 于 上 缘 板 沿 水 平 向 周围 突出 , 糊 
状 区 横向 周边 的 上 方 的 冷却 速率 突然 减 小 , 铸件 
表面 凝固 前 沿 推移 减 慢 , 因而 顶端 的 下 四 程度 减 
小 . 综 上 所 述 , 使 得 KK 曲线 呈 S 形 . 引 唱 1,2 和 3 
的 Ks 依次 略 有 减 小 . 这 是 由 于 叶 身 与 引 唱 块 连接 
方式 和 激 冷 面 接触 面积 的 不 同 , 对 铸件 心 部 的 许 
固 前 沿 推移 速率 的 影响 强 弱 也 就 不 同 . 引 唱 1 2 
和 3 的 上 凸 阶 段 较 短 , 下 凹 阶段 的 下 止 程度 很 大 ， 
趋势 上 容易 导致 晶 粒 反 发 散 生长 . 引 唱 4 与 激 冷 盘 
的 接触 面积 大 且 直 接连 接 叶 身 , 因此 引 晶 4 的 K. 明 
显 小 于 前 3 种 引 晶 . 
4.1.2 不 同 引 晶 方 式 的 唱 粒 组 织 分 析 图 9 
为 引 晶 1 下 拉 速 7 mm/min 时 计算 和 实验 得 到 的 叶 
身 晶 粒 组 织 . 模拟 和 实验 结果 均 显 示 : 引 晶 块 内 的 
鲁 粒 经 竞争 生长 后 得 到 大 致 平行 于 竖 直 方向 的 柱 
状 唱 , 然后 进入 下 缘 板 横向 生长 , 最 后 只 有 少 部 分 
进入 叶 身 , 形成 过 宽 的 晶 粒 , 叶 贫 可 见 2 个 宽度 超 
过 弦 线 30% 的 晶 粒 (如 图 9 中 I 所 示 ), 从 进 气 边 到 叶 
盆 侧 30 mm 宽度 和 到 叶 背 侧 45 mm 宽度 的 范围 几 
乎 都 被 同一 个 晶 粒 占据 . 排 气 边 背 侧 离 引 晶 块 相对 
较 近 , 因而 靠 排 气 边 唱 粒 数 相 对 多 一 些 , 唱 粒 宽度 
比 进 气 边 略 小 (如 图 9 中 开 所 示 ). 叶 盆 下 段 靠 进 气 
边 的 部 分 , 出 现 了 明显 的 带 有 横向 唱 界 的 断 唱 (如 
图 9 中 II 所 示 ). 晶 粒 总 体 上 略 有 反 发 散 趋势 , 使 得 
叶 背 出 现 大 块 三 角形 露头 唱 ( 如 图 9 中 IV 所 示 ). 
图 10 为 引 晶 2 下 拉 速 7 mmmin 时 计算 和 实验 
得 到 的 叶 身 唱 粒 组 织 , 比 引 晶 1 有 所 改善 . 模拟 和 实 
验 结 果 均 显示 : 增加 引 唱 块 后 , 进入 叶 盆 的 唱 粒 数 
增加 , 唱 粒 宽度 减 小 , 但 靠 进 气 边 部 分 的 晶 粒 仍 比 
靠 排 气 边 部 分 宽 . 盆 、 背 两 面 靠 排 气 边 出 现 少许 细 
窄 断 晶 (如 图 10 中 I 所 示 ), 叶 背 靠 进 气 边 出 现 大 块 
三 角形 晶 粒 (如 图 10 中 开 所 示 ), 但 与 上 下 缘 板 相连 ， 
不 是 断 晶 . 叶 贫 下 段 靠 进 气 边 的 带 有 横向 晶 界 的 断 
晶 大 大 减少 (如 图 10 中 II 所 示 ). 晶 粒 总 体 上 略 有 反 
发 散 趋势 . 
图 11 为 引 晶 3 下 拉 速 7 mm/min 时 计算 和 实验 
得 到 的 叶 身 晶 粒 组 织 , 比 引 晶 2 有 所 改善 . 模拟 和 实 


才 
\ 


呈 呈 所 


202303.00464v1 


chinaXiv 


506 金 属 


ChinaXiv 合 作 期 刊 


学 报 第 51 郑 


9 引 晶 方式 1 下 拉 速 7 mm/min 时 得 到 的 叶 盆 和 叶 背 晶 粒 组 织 的 模拟 结果 和 实验 结果 
Fig.9 Numerical (a, c) and experimental (b, d) morphologies of blade basin (a, b) and blade back (c, d) obtained by using 


starter 1 at withdrawal rate 7 mmy/min ( I—excessively wide grains; II—narrower grains; II 一 fragmented grains; 


IV 一 triangular grains) 


10 引 唱 方式 2 下 拉 速 7 mm/min 时 得 到 的 叶 贫 和 叶 背 唱 粒 组 织 的 模拟 结果 和 实验 结 
Fig.10 Numerical (a, c) and experimental (b, d) morphologies of blade basin (a, b) and blade back (c, d) obtained by using 


starter 2 at withdrawal rate 7 mm/min (I—fragmented grains; I[—triangular grains; II 一 全 agmented grains) 


验 结果 均 显 示 : 进入 叶 身 的 唱 粒 进一步 增多 , 这 是 
由 于 从 引 蝇 段 生 长 的 蝇 粒 可 以 直接 长 到 叶 身 , 而 不 
是 通过 缘 板 , 这 样 避免 了 缘 板 处 晶 粒 生长 到 叶 身 时 
晶 粒 取向 的 偏离 和 晶 粒 粗大 . 叶 盆 进 气 边 侧 与 排 气 


边 侧 的 晶 粒 宽度 相差 很 小 , 叶 背 不 再 出 现 大 三 角 唱 
粒 , 而 叶 背 中 部 出 现 1~2 个 分 又 状 窗 唱 粒 (如 图 11 中 
I 所 示 ), 其 两 侧 各 有 一 个 略 宽 的 晶 粒 . 叶 盆 下 段 靠 两 
边 的 部 分 都 有 多 股 分 又 的 树 状 晶 粒 (如 图 11 中 开 所 
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示 ). 唱 粒 总 体 上 有 反 发 散 趋势 . 分 又 唱 粒 会 对 表面 
横 截 线 唱 粒 数 的 统计 造成 干扰 , 因为 同一 个 品 粒 对 
应 很 多 个 线段 , 不 能 真实 反映 品 粒 生长 情况 . 因此 ， 
生产 中 应 尽量 减少 分 又 晶 粒 . 

原 有 3 种 引 晶 方式 下 计算 的 典型 晶 粒 的 形状 和 
位 置 都 与 实验 吻合 ,证 明了 计算 的 可 靠 性 . 

12 为 引 晶 4 下 拉 速 7 mm/min 时 计算 得 到 的 
叶 身 晶 粒 组 织 . 进入 叶 盆 和 叶 背 的 唱 粒 数 都 有 大 幅 
度 增加 , 品 粒 宽度 大 大 减 小 , 叶 倪 、 叶 背 进 气 边 侧 与 
排 气 边 侧 的 唱 粒 宽度 都 相差 很 小 , 叶 贫 、 叶 背 下 段 
都 出 现 少量 罕 细 露头 唱 , 但 没有 明显 的 横向 唱 界 ， 
叶 背 靠 排 气 边 部 分 有 个 别 三 角 型 露头 晶 粒 . 没有 出 
现 分 又 状 晶 粒 . 由 于 糊 状 区 下 四 程度 减 小 , 晶 粒 总 


@) 


体 上 沿 竖 直方 向 平行 , 基本 没有 发 散 趋势. 

不 同 引 蝇 方 式 得 到 的 叶 身 Kxo 的 变化 如 图 13 所 
示 . 可 见 , 唱 粒 数 随 着 凝固 前 沿 推移 而 减少 , 引 品 1， 
2 和 3 的 Ko 依次 增加 . 引 唱 4 的 唱 粒 数 明显 高 于 前 3 
种 引 唱 . 引 蝇 4 叶 身 下 端 Kxo 下 降 较 快 , 表示 断 唱 较 
多 , 但 如 前 文 所 述 , 都 是 窗 小 露头 晶 , 横向 晶 界 较 少 . 

不 同 引 蝇 方 式 得 到 的 叶 身 Kswo 变 化 如 图 14 所 
示 . 可 见 , 引 唱 1,2 和 3 的 Kso 依 次 增加 , 引 唱 1 和 2 
的 Ksxo 随 凝固 前 沿 推 移 的 变化 规律 与 Kvo 一 致 , 而 引 
品 3 的 Kswo 在 叶 身 下 半 段 , 先 下 降 后 回升 , 围绕 一 个 
平台 值 上 下 波动 , 进入 叶 身 上 半 段 后 继续 下 降 . 这 
主要 来 源 于 叶 身 下 半 段 的 树 状 品 粒 的 干扰 , 同一 个 
树 状 晶 粒 在 Ksxo 的 计算 中 对 应 很 多 个 线段 , 使 得 


11 引 唱 方式 3 下 拉 速 7 mm/min 时 得 到 的 叶 盆 和 叶 背 晶 粒 组 织 的 模拟 结果 和 实验 结果 
Fig.11 Numerical (a, c) and experimental (b, d) morphologies of blade basin (a, b) and blade back (c, d) obtained by using 


starter 3 at withdrawal rate 7 mm/min (1, I[—branching grains) 


I 


(b) 


12 引 晶 方式 4 下 拉 速 7 mm/min 时 得 到 的 叶 贫 和 叶 背 唱 粒 组 织 的 模拟 结果 
Fig.12 Numerical morphologies of blade basin (a) and blade back (b) obtained by using starter 4 at withdrawal rate 


7 mm/min 
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Sz 


Kswo 升 高 . 引 晶 4 


出 现 分 又 晶 粒 引 起 的 Kswo 回 升 . 


证 了 模型 和 参数 的 正确 


的 Kswo 高 于 前 3 种 引 晶 , 并 且 没 有 


原 有 的 3 种 工艺 的 模拟 结果 都 与 实验 吻合 , 验 
性 . 由 模拟 结果 可 知 , 增 大 


铸件 与 激 冷 面 的 接触 面积 有 利于 减 小 糊 状 区 宽度 
和 下 凹 度 、 增 加 进入 截面 的 晶 粒 数 、 减 少 断 唱 和 过 


宽 的 品 粒 . 


大 . 引 晶 1,2 和 3 的 到。 整体 都 1 


本 节 研 究 发 现 ,在 4 种 引 唱 方式 中 , 引 晶 3 和 4 
较 优 , 引 唱 3 尚 存 在 唱 粒 反 发 散 问题 , 这 可 以 通过 优 
化 拉 速 来 解决 . 

4.2 基于 数值 模拟 的 抽 拉 速度 优化 
如 前 文 所 述 , 引 唱 1, 2 和 3 中 , 同一 拉 速 下 , 引 
晶 方 式 对 K 有 和 Ks 影响 不 大 , 而 对 唱 粒 数 的 影响 很 


i 大 , 品 粒 都 有 反 发 


散 趋势 , 应 考虑 减 小 拉 速 . 生产 经 验 表 明 , 其 他 参数 


固 


时 , 才 会 对 品 粒 有 显著 的 影响 . 因此 , 本 工作 采 


定 的 情况 下 , 不 同 工 艺 的 拉 速 相差 2 mm/min 以 上 
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对 于 引 晶 3, 计算 了 3,5 和 7 mmmin 下 的 温度 
场 和 微观 组 织 , 为 验证 计算 的 准确 性 , 对 引 晶 3 在 拉 
速 为 5mm/min 下 进行 了 实验 . 

4.2.1 不 同 拉 速 下 的 糊 状 区 形态 分 析 图 15 
为 引 晶 3 和 4 在 不 同 拉 速 得 到 的 叶 身 KK 变化 . 可 见 ， 
在 叶 身 下 半 段 , 拉 速 对 KK 的 影响 很 小 , 同一 引 晶 方 
式 的 Ku 几乎 重合 . 这 是 由 于 凝固 前 沿 离 激 冷 面 较 近 
时 , 激 冷 盘 的 散热 效率 决定 了 糊 状 区 中 的 温度 梯 
度 , 而 且 模 壳 温 度 受 激 冷 盘 的 影响 很 大 , 挡 板 相对 
位 置 对 模 壳 温度 的 影响 被 弱化 , 因此 , 抽 拉 过 程 中 
挡 板 相 对 位 置 的 变化 对 糊 状 区 温度 梯度 的 影响 就 
更 小 . 此 时 , K, 主 要 取决 于 不 同 引 晶 方 式 下 铸件 与 
激 冷 盘 的 接触 面积 . 在 叶 身 下 半 段 , K, 随 着 拉 速 的 
增加 而 有 所 增 大 . 
图 16 为 引 唱 3 和 4 在 不 同 拉 速 得 到 的 叶 身 及。 
变化 . 可 见 , 叶 身 下 端的 Kw 与 工艺 参数 的 关系 类 似 


不 同 拉 速 之 间 至 少 相 差 2 mm/min. 


80 一 
| 一 一 Starter1,7 Mmm/min 
70 上 上 ---- Starter 2, 7 mm/min 
二 Starter 3, 7 mm/min 
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50| ~ 
© ， h 
Xe 40F. 


30 上 ~ 
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13 不 同 引 蝇 方 式 下 叶 身 横 截面 晶 粒 数 Ko 沿 高 度 的 
变化 
Fig.13 Variations of grain number on lateral section Kxo in 


blade body with height for different starters 


80 一 
一 一 Starter1,7 mm/min 
70 ia ---- Starter 2, 7 mm/min 
和 Starter 3, 7 mM/min 
60 上 一 Starter 4, 7 mm/min 
50 上 上 1 { Py 


Height of solidification front / mm 
14 不 同 引 晶 方 式 下 叶 身 横 截 线 晶 粒 数 Kswo 沿 高 度 的 
变化 
Fig.14 Variations of surface grain number on lateral line 


Ksvo in blade body with height for different starters 


的 


于 Ks, 受 拉 速 影 响 很 小 . 从 距 下 缘 板 约 12 mm 处 开 
始 , Kew, 随 着 拉 速 增加 而 显著 增加 , 且 拉 速 的 对 Ko 
的 影响 不 断 增 大 . 对 于 同一 种 引 蝇 , 拉 速 降低 时 , 其 
总 体 下 四 程度 降低 , 上 凸 阶 段 延 长 ,下 四 阶段 变 短 . 


10 
一 一 Starter3, 3 mm/min 
9 Le Starter 3, 5 mm/min 
8 Starter 3, 7 mM/min 
7L 一 Starter 4, 5 mm/min 
一 Starter 4, 7 mm/min 
E 6 : 
四 
~ 5F 
Sl 
上 
2 上 
1 上 
0 1 | | 
40 50 60 70 80 90 100 


Height of solidification front / mm 


15 引 晶 3 和 引 晶 4 不 同 拉 速 下 叶 身 KK 沿 高 度 的 变化 
Fig.1S Variations of Ku in blade body with height at differ- 


ent withdrawal rates for starters 3 and 4 


0.4 
一 一 Starter3, 3 mm/min 
03| ---- Starter 3, 5 mm/min 
ee Starter 3, 7 mm/min 
ie Starter 4, 5 mm/min 
“| --- Starter 4, 7 mm/min 
5 0.1 
x 
0.0 上 一 
-0.1 
:2 二 
1 1 上 


40 50 60 70 80 90 100 
Height of solidification front /mm 
16 引 虽 3 和 引 唱 4 不 同 拉 速 下 叶 身 Ke 沿 高 度 的 变化 
Fig.16 Variations of K.,, in blade body with height at dif- 


ferent withdrawal rates for starters 3 and 4 
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4.2.2 不 同 拉 速 下 的 晶 粒 组 织 分 析 17 为 
引 晶 3 下 拉 速 5 mm/min 时 计算 和 实验 得 到 的 叶 身 
晶 粒 组 织 . 模拟 和 实验 结果 均 显示 : 进入 叶 身 的 晶 
粒 数 、 晶 粒 宽度 都 与 拉 速 为 7mnymin 时 相差 不 大 ， 
叶 背 靠 进 气 边 的 部 分 出 现 较 宽 的 晶 粒 (如 图 17 中 I 
所 示 ) 和 三 角 晶 粒 (如 图 17 中 开 所 示 ), 但 宽度 没有 超 
过 30% 弱 长 , 叶 背 中 间 和 叶 盆 人 靠 进 气 边 部 分 的 分 又 
唱 粒 减少 . 最 重要 的 是 , 唱 粒 总 体 反 发 散 趋势 比 拉 
速 为 7 mnymin 时 减 小 , 特别 是 改善 了 叶 贫 部 分 的 唱 
粒 平行 度 . 

18 为 引 唱 3 下 拉 速 3 mm/min 时 计算 得 到 的 
叶 身 晶 粒 组 织 . 模拟 结果 显示 : 进入 叶 身 的 晶 粒 
数 、 唱 粒 宽度 都 与 拉 速 为 7 mm/min 时 相差 不 大 , 叶 
背 出 现 倒 三 角 唱 粒 , 叶 盆 靠 进 气 边 的 部 分 有 多 股 分 


又 的 树 状 晶 粒 . 晶 粒 总 体 上 时 发 散 趋势 . 

19 为 引 唱 3 和 4 在 不 同 拉 速 下 得 到 的 叶 身 
Kwo 变 化 . 可 见 , 品 粒 数 随 着 凝固 前 沿 推移 而 减少 , 
拉 速 对 Kro 的 影响 较 小 . 

图 20 为 引 晶 3 和 4 在 不 同 拉 速 下 得 到 的 叶 身 
Kswo 变 化 . 可 见 , Kswo 随 着 拉 速 的 增加 而 增加 . 引 品 3 
在 不 同 拉 速 下 , Kswo 都 在 叶 身 下 段 距 下 缘 板 约 10 mm 
处 明显 回升 . 

Ksvo 回 升 说 明 引 晶 3 很 容易 形成 树 状 唱 粒 , 其 
形成 趋势 受 拉 速 的 影响 不 大 . 对 品 粒 进行 三 维 解剖 
分 析 , 发 现 是 由 于 叶 身 的 晶 粒 有 2 个 来 源 : 从 下 缘 板 
凸 棱 处 引 唱 块 进入 下 缘 板 , 然后 经 过 横向 竞争 生长 
进入 叶 身 的 唱 粒 ; 从 叶 盆 下 方 引 晶 块 经 过 纵向 竞争 
生长 直接 进入 叶 身 的 唱 粒 . 横向 和 纵向 两 组 唱 粒 交 


17 引 晶 方式 3 下 拉 速 5 mm/min 时 得 到 的 叶 盆 和 叶 背 唱 粒 组 织 的 模拟 结果 和 实验 结 
Fig.17 Numerical (a, c) and experimental (b, d) morphologies of blade basin (a, b) and blade back (c, d) obtained by using 


starter 3 at withdrawal rate 5 mm/min (I—wide grains; II—triangular grains) 


18 引 晶 方式 3 下拉 速 3 mm/min 时 得 到 的 叶 贫 和 叶 背 唱 粒 组 织 的 模拟 结果 
Fig.18 Numerical morphologies of blade basin (a) and blade back (b) obtained by using starter 3 at withdrawal rate 
3 mm/min 
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19 引 晶 3 和 引 晶 4 不 同 拉 速 下 叶 身 Kxo 沿 高 度 的 变化 

Fig.19 Variations of Kxo in blade body with height at dif- 


ferent withdrawal rates for starters 3 and 4 
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20 引 晶 3 和 引 晶 4 不 同 拉 速 下 叶 身 Kswo 沿 高 度 的 变化 
Fig.20 Variations of Ksvo in blade body with height at dif- 


ferent withdrawal rates for starters 3 and 4 


汇 进 入 叶 身 时 都 大 致 沿 择优 取向 , 相互 争夺 叶 身 的 
生长 空间 , 再 加 上 叶 身 复杂 曲面 的 干扰 , 有 时 会 使 
得 2 个 唱 粒 交 蔡 占据 空间 , 从 而 使 其 中 一 个 唱 粒 被 
另 一 个 晶 粒 割裂 为 多 个 分 支 , 形成 树 状 唱 粒 . 

引 品 1 只 有 下 缘 板 凸 校 处 引 晶 ; 引 唱 2 添 加 的 
引 章 块 中 的 唱 粒 需 先进 入 下 缘 板 , 与 来 自 下 缘 板 
凸 楼 处 的 晶 粒 一 起 在 缘 板 中 进行 横向 竞争 生长 ， 
优势 晶 粒 进入 叶 喘 ; 引 品 4 添加 的 全 弦 线 引 品 块 
中 经 过 纵向 竞争 生长 的 唱 粒 进入 叶 身 时 占据 主 
导 , 缘 板 横向 生长 的 唱 粒 难以 进入 叶 身 . 总 之 , 引 
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是 5 mm/min, 唯一 的 缺点 是 分 又 晶 粒 难以 避免 . 

拉 速 7mmmmin 下 引 唱 4 的 唱 粒 数 和 唱 粒 宽度 
优 于 前 3 种 引 晶 , 而 和 Ks 与 前 3 种 引 蝇 明显 不 
同 , 有 必要 再 研究 引 晶 4 下 不 同 拉 速 的 影响 , 特别 是 
5 mm/min 对 引 唱 4 来 说 是 否 为 最 佳 拉 速 . 此 外 , 3 
鲁 4 的 唱 粒 没有 反 发 散 趋势 , 其 Ru<0 的 阶段 占 叶 
身高 度 的 近 2/3, 因此 , 引 晶 4 再 提高 拉 速 对 优化 
晶 粒 组 织 作 用 不 大 . 因此 , 将 引 晶 4 的 拉 速 减 小 为 
5 mm/min, 进行 了 相应 工艺 的 宏 微观 模拟 计算 , 并 
与 拉 速 7 mnymin 时 进行 比较 . 

由 图 15 可 见 , 引 唱 4 不 同 拉 速 的 影响 与 引 晶 3 
相似 , 但 引 晶 方式 的 影响 比 拉 速 的 影响 强 , 引 晶 4 最 
大 拉 速 下 的 K& 仍 小 于 引 晶 3 最 小 拉 速 下 的 K.. 

由 图 16 可 见 , 引 唱 4 不 同 拉 速 的 影响 与 引 唱 3 

相似 . 由 于 激 冷 面 接触 面积 的 影响 , 引 晶 4 的 下 四 程 
度 显 著 低 于 引 晶 3. 引 晶 4 在 拉 速 为 7 mm/min 的 KK 
曲线 与 引 晶 3 拉 速 为 mm/min 的 相 比 , 比较 接近 ， 
但 是 , 前 者 的 绝对 值 在 大 部 分 阶段 比 后 者 小 , 亦 即 
前 者 在 上 凸 阶段 的 上 凸 程度 比 后 者 小 , 下 四 阶段 的 
下 四 程度 也 小 于 后 者 , 前 者 更 接近 平面 且 变 化 幅度 
更 小 . 由 前 文 可 知 在 引 晶 3 下, 5 mm/min 是 较 优 拉 
速 . 因此 , 引 晶 4 下 7 mm/min 有 可 能 是 最 佳 拉 速 . 
图 21 为 引 晶 4 下 5 mmymin 时 计算 得 到 的 叶 身 
晶 粒 组 织 . 模拟 结果 显示 , 其 叶 贫 的 晶 粒 数量 比 拉 
速 为 7 mmmmin 时 少 , 叶 身 下 段 的 露头 晶 比 拉 速 为 
7 mm/min 时 少 , 但 进 气 边 上 有 一 个 断 晶 . 叶 背 中 间 
部 分 出 现 2 个 分 又 晶 粒 . 晶 粒 总 体 略 有 发 散 趋 势 , 特 
别 是 叶 背 靠 排 气 边 的 晶 粒 , 严重 倾斜 . 因此 , 引 晶 4 
下 拉 速 为 5 mm/min 的 唱 粒 组 织 比 7 mm/min 时 差 . 

由 图 19 和 20 可 见 , 引 晶 4 在 2 种 拉 速 下 得 到 的 
叶 身 Ko 和 Kso 都 明显 大 于 引 晶 3, 且 拉 速 为 7 mm/ 
min 时 晶 粒 数 大 于 5 mm/min. 

引 晶 4 的 温度 场 和 晶 粒 模拟 表明 其 较 优 拉 速 为 
7 mm/min, 而 且 其 各 项 评价 指标 均 优 于 或 接近 于 引 
晶 3 在 最 优 拉 速 5mm/min 下 的 结果 , 并 且 避 免 了 引 


一 二 


唱 1,2 和 4 都 难以 在 叶 身 下 段 出 现 横向 和 纵向 两 
组 晶 粒 竞争 生长 的 情况 , 所 以 不 容易 形成 树 状 唱 
粒 . 在 实验 结果 中 , 也 只 有 引 晶 3 的 叶 身 观察 到 了 
分 又 状 的 晶 粒 . 

前 述 研究 发 现 , 拉 速 对 和 Ks 影响 很 大 , 而 对 
截面 晶 粒 数 和 晶 粒 宽度 的 影响 很 小 . 引 唱 1 和 2 的 
唱 粒 数 和 晶 粒 宽度 都 不 佳 , 因此 不 再 进行 引 晶 1 和 2 
的 拉 速 优化 . 

引 晶 3 的 温度 场 和 唱 粒 模拟 表明 , 其 最 佳 拉 速 


由 


唱 3 的 分 又 晶 粒 . 
4.3 不 同 工 艺 的 综合 比较 

综合 分 析 所 有 的 计算 结果 , 引 唱 方式 1 和 2 在 
拉 速 为 7mm/min 下 的 糊 状 区 宽度 . 唱 粒 宽度 、 唱 粒 
数目 不 合格 , 断 唱 产生 的 横向 唱 界 较 明 显 . 而 前 述 
研究 表明 , 多 数 引 晶 方式 都 可 以 通过 改变 拉 速 得 到 
较 好 的 晶 粒 平行 度 , 但 拉 速 对 其 他 指标 的 改善 不 显 
著 . 因此 , 引 唱 方式 1 和 2 在 其 他 拉 速 下 也 不 会 得 到 
优化 的 晶 粒 组 织 . 
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(b) 


21 引 品 方式 4 下 拉 速 S mm/min 时 得 到 的 叶 贫 和 叶 背 唱 粒 组 织 的 模拟 结果 
Fig.21 Numerical morphologies of blade basin (a) and blade back (b) obtained by using starter 4 at withdrawal rate 


5 mm/min 


而 对 引 晶 方式 3 和 4 的 研究 表明 , 它们 最 优 拉 
速 分 别 是 5 和 7 mm/min, 在 优化 拉 速 下 都 得 到 了 比 
引 晶 方 式 1 和 2 更 优 的 组 织 . 

可 见 , 引 唱 方式 1 和 2 应 该 据 弃 . 如 果 后 续 高 
精度 切削 的 成 本 不 高 , 可 选择 引 唱 方式 4, 拉 速 为 
7 mm/min, 以 获得 最 优 的 微观 组 织 和 最 高 的 铸造 效 
率 ; 也 可 选择 引 晶 方式 3, 拉 速 为 5 mm/min, 以 获得 
相对 较 好 的 微观 组 织 . 

5 结论 

(1) 建立 了 宽 弦 航空 叶片 定向 凝固 宏 、 微 观 多 尺 
度数 学 模型 , 基于 自主 开发 的 程序 对 不 同 工 艺 的 温 
度 场 和 唱 粒 组 织 进行 了 模拟 , 计算 结果 与 实验 结 
吻合 较 好 , 铸件 表面 大 部 分 晶 粒 的 位 置 和 形状 一 一 
对 应 . 

(2) 建立 了 宏 微 观 数学 评价 判 据 , 实现 了 不 同 工 
艺 下 糊 状 区 宽度 、 凝 固 前 沿 凸 止 度 、 唱 粒 数 的 定量 
比较 , 获得 了 工艺 参数 对 糊 状 区 和 晶 粒 组 织 的 影响 
规律 . 

(3) 引 唱 块 与 激 冷 盘 接 触 面 积 变化 不 大 时 , 引 
唱 方 式 主要 影响 晶 粒 数 和 晶 粒 宽度 , 拉 速 主要 影响 
糊 状 区 宽度 和 凝固 前 沿 凸 凹 度 ; 拉 速 增 大 时 下 四 
趋势 越 大 .下 止 阶段 变 长 , 糊 状 区 宽度 略 有 增加 , 反 

(4) 铸件 与 激 冷 盘 之 间 的 导热 面积 增 大 时 , 糊 状 
区 变 窗 , 上 凸 趋势 变 大 . 凝固 前 沿 离 激 冷 面 较 近 时 ， 
对 该 面积 的 影响 较 明 显 , 拉 速 的 影响 被 弱化 . 

(5) 基于 数值 模拟 和 判 据 , 研究 了 不 同 引 唱 方式 
和 拉 速 下 晶 粒 组 织 的 变化 机 理 , 发 现 原 有 的 3 种 引 
唱 方 式 中 的 最 优 者 为 叶 盆 下 方 引 晶 , 并 确定 了 相应 
的 最 佳 拉 速 . 

(6) 基于 模拟 预测 结果 和 不 同 工 艺 影响 下 的 微 
观 组 织 形 成 理论 , 提出 了 新 的 全 弦 线 引 晶 方法 , 3 


对 其 拉 速 进行 了 优化 , 可 获得 枝 晶 更 细密 、 唱 粒 更 
平行 .柱状 唱 更 多 的 铸件 , 使 生产 效率 进一步 提高 . 
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